DVH computation comparisons using conventional methods and optimized FFT algorithms for brachytherapy

In anatomy based optimization procedures for large volume implants the calculation of dose-volume histograms (DVH) accounts for the major part of the time involved and can be as long as a few hours. This time is proportional to the number of seeds or source dwell positions required for the implant. A procedure for the calculation of brachytherapy seed dose distribution calculation employing fast Fourier transforms (FFT) and the convolution theorem has been described by others and was supposed to significantly improve the speed of the dose distribution computation. Using new significantly improved FFT algorithms and various other optimization techniques we have compared the calculated differential and integral DVHs in HDR brachytherapy with a single stepping source using actual clinical implants. This is so that we could assess the efficiency and accuracy of the FFT method with that of conventional methods. Our results showed that the FFT based method of calculating DVHs in brachytherapy is comparable in speed with conventional dose calculation methods, but only for implants with more than 287 sources. It is therefore of limited practical use even for large implants. This result is in direct opposition to the claim by other authors.